Estimation of gravitational bending of light in the weak deflection limit 
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After discussing some subtle issues concerning the computation of deflection angle, a general but 
simple expression of bending angle of light rays in weak defelction limit has been presented for 
a general static and spherically symmetric space-time. In this context the importance of proper 
choice of polar axis has been highlighted. Applying the prescribed method the explicit expression of 
bending angle up to an accuracy of second order in mass is obtained for the Schwarzschild geometry 
without restricting the source and the observer to be at infinite distance away from the lens. 
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I. INTRODUCTION 

Light rays while passing a massive body are bent to- 
wards that mass. The prediction of general relativity in 
respect to the angle of deflection for the light ray just 
grazing the limb of the Sun was confirmed in 1919 when 
the apparent angular shift of stars close to the limb of the 
sun was measured during a total solar eclipse. This grav- 
itational deflection angle has since been repeatedly con- 
firmed with high-precision measurements and it provides 
compelling evidence in favor of general relativity. An im- 
portant application of the deflection effect concerns the 
so-called gravitational lens that acts as a tool to study 
various features of the universe including estimation of 
masses of galaxies [1] . 

Conventionally the quantum of bending of light rays 
due to a spherically symmetric mass distribution (lens) 
is computed by studying the trajectory of light rays in the 
gravitational field of the lens. Usually light trajectory is 
obtained from the null geodesic equation either using per- 
turbative approach (see for instance [2-3]) or by directly 
integrating the null geodesic equation for the orbit equa- 
tion of light [4-5] . When both the source and the observer 
are situated far far away from the lens in compare to the 
distance of closest approach of light rays from the lens, 
the angle between the two asymptotes of the light orbit is 
considered as the bending angle though such an equality 
is not apparent from the lensing diagram. More impor- 
tantly the justification for such an equivalence has not 
commonly been discussed in the standard monographs 
of general relativity. As a result several issues concern- 
ing estimation of bending angle remain unclear to many 
and the equality of bending angle and the angle between 
the two asymptotes of the light trajectory has often been 
used in situations where it is not truly valid. The ex- 
istence of different methods of estimating bending an- 
gle has further enhanced the perplexity of the problem. 
Consequently different expressions of bending angle for 
the same lensing problem, particularly whenever gravita- 
tional bending is considered in situations other than the 
classic case of infmctly distant observer and source, have 
been found in the literature that creats greater confusion. 

The purpose of the present article is to address such 



subtle issues concerning computation of bending angle. 
We will particularly highlight the importance of the 
proper choice of polar axis for framing the deflection sce- 
nario and calculating the deflection angle. Subsequently 
we will describe a straightforward and general approach 
of computing bending angle when the source and /or ob- 
server are located at a finite distance away from the lens. 
After providing a rigorous analytical framework of esti- 
mating bending angle for the general static spherically 
symmetric metric we will calculate the explicit expres- 
sions of bending angle for the Schwarzschild geometry. 

The organization of the article is follows. In the next 
section we would discuss the null geodesic equations for 
the general static spherically symmetric metric. In sec- 
tion 3 we would write down the orbit equation of light in 
the Schwarzschild geometry. The deflection angle will be 
computed in section 4 and finally we would conclude in 
section 5. 



II. EQUATION OF MOTION FOR PHOTONS 

We consider the general static and spherically symmet- 
ric spacetime which is given by (we use geometrized units 
i.e. G = l, c=l) 



ds 2 = -B(f)dt 2 + A(f)df 2 + r 2 D{f) (d9 2 + sin 2 6d(j) 2 ) 

(1) 
In principle one may recast the above metric to 



the standard form ds 2 = —f(r)dt 2 



g(r)dr 



(d9 2 + sin 2 9d(f) 2 ) just by transforming to a new ra- 
dial variable r where r = fy/D(f), but in some cases no 
explicit expression for the transformation equation exists 
or the transformation equation is too complex to deal 
with. Hence at this stage we would continue with the 
form as given by Eq. (1). 

Instead of obtaining the geodesic equations directly 

,2 A 

from the standard equations for a geodesic ^^r + 

^u.vlt~lt~ ~ 0' one ma y § et tnose from the Euler La- 
grange kind of equations, viz. [2] 



dp \ dx^ J dx^ 



(2) 



where Q = g^ v dx IJ 'dx v is the square of the Lagrangian 
(L) of the gravitational field, dx^ = dx^/dp, p is some 
affinc parameter along the geodesic x^(p). 

The square of the Lagrangian for the general metric 
(1) is given by 



C = -B(f)i 2 + A(f)f 2 + f 2 D(f) (e 2 + sin 2 9d^ 2 \ (3) 
The resulting four geodesic equations are given by 



dp V dp) 



(4) 




dp \ dp 
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(5) 



(6) 



where E is a constant of motion. Before proceeding 
further we would see what these constant of motion re- 
fer to. The Eqs. (1), (8), (9) and (10) together give 
ds = dp/ E. For photon obviously 1/E = whereas for 
massive particle E > 0. Since £ does not depend explic- 
itly on t and <f>, the t and <j> components (pt and p^ ) of 
a particle's 4-momentum (p^ oc dx^ /ds) are conserved 
along geodesies. For an observer at rest at infinity the 
energy per unit rest mass of a particle is e = p t . There- 
fore, we have from Eq.(4) e = g tt p t = E. On the other 
hand p,p = g^d^/ds — f 2 Dd(f>/ ds = bE. Thus b repre- 
sents the impact parameter at large distances. Here note 
that the unit of the constants of motion depends on the 
choice of afine parameter and the normalization of the 
time coordinate. For instance if one takes dp = ds, then 
b will be angular momentum per unit mass (or for mass 
less particle it is the angular momentum), E = 1 and k 
(of Eq.(8)) would then describe the energy per unit rest 
mass. 

From Eqs. (8), (9) and (10), we finally obtain for pho- 
tons 



A 



dr 



D 
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= 



(11) 



Integrating the above equation we will get as a func- 
tion off for the light orbit. However, to describe an orbit 
one usually expresses radial variable as a function of an- 
gular coordinate. Differentiating Eq.(ll) with respect to 
4> one gets 



± (r 2 D d 4- 1 = o 

dp \ dp, 



(7) 



where primes denotes differentiation with respect to f . 
If we choose 9 = tt/2 and d6/dp = initially, Eq. (6) war- 
rants that they would remain the same always. Thus for 
simplicity hereafter we would consider the motion only 
in the equatorial plan. 

The Eq. (4) straightway gives 
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(12) 



where u = 1/f. For a given metric i.e. for given B, A 
and D the solution of the above second order differential 
equation gives the photon orbit. 



III. THE ORBIT EQUATION FOR PHOTON IN 
THE SCHWARZSCHILD GEOMETRY 



dt 
dp 



k_ 
B 



(8) 



where fc is a constant. Normalizing time coordinate 
suitably we can set k = 1. Integrating Eq.(7) 



The first step of calculating bending angle is to obtain 
the orbit equation. We would discuss hereunder both the 
popular approaches of solving geodesic equations for the 
light trajectory. 



f 2 D-f- = b 
dp 



(9) 



where b is a constant of integration. The Eq.(5) to- 
gether with Eqs. (8), (9) give 



A. Perturbation method 

Let us consider first the case when there is no lens and 
hence the geometry is flat (Euclidean) i.e. A = B = D = 
1. In such a situation the Eq.(12) reduces to 



dp) ~B + f 2 ~D~E 2 ~ 



(10) 



d 2 u 



(13) 



Hence the orbit equation of undeficcted (straight line) 
light ray in flat space not containing the origin of the 
coordinate system is given by 



u= —sin(4> + t;) 



(14) 



where R is the perpendicular distance from the origin 
of the co-ordinate system to the path of the light rays 
and £ is the angle that the light rays made with the polar 
axis at the point of intersection. The angle <f> is measured 
from the centre of the co-ordinate system with respect to 
the polar axis. Note that in the absence of any lens one 
is free to take any point as the origin of the coordinate 
system and also may set any direction of his/her choice 
as the direction of the polar axis. 

Now we consider the propagation of light ray in the 
gravitational field of the lens that describes a spherical 
distribution of mass. We are not considering the rota- 
tional motion of the lens/observer. Since we have to deal 
with the gravitational field of the Lens it is much conve- 
nient to take the Center of the Lens (L) as the origin of 
the coordinate system. 

To have a solution of the Eq.(12) one has to supply 
explicit expressions for A, B and D. Now referring to 
the case of the Schwarzschild geometry 
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(18) 



The above equation represents the orbit equation for 
the light ray up to an accuracy of m 2 . Since the shape 
of the orbit has changed due to the presence of the lens 
the distance of closest approach is also likely to be al- 
tered. Differentiating both sides of the above equation 
with respect to 4> one gets 



1 dr 



cos((f> + £) 
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15m 2 , 
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(19) 



At the point of closest approach & vanishes which 
occurs at <j> + £ = n/2. Thus the distance of closest 
approach r is given by 



2,„ -Vl^.J^ 



(de 2 
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ds 2 = -(1 )dt 2 +{\ 

r r 

(15) 

i.e. when B = 1 — 2m/r, A = 1/5, D = 1 and r = r 

the null geodesic equation (12) reduces to 



d 2 u 



3mu 



(16) 



When the gravitational field is weak we can expect that 
the light path will deviate by a small amount from the 
undeficcted (straight line) path and accordingly in such 
a situation we will look for a perturbative solution about 
the undeficcted orbit. Therefore, we may take the solu- 
tion as u = u +u\, where u is the undeflectcd trajectory 
as given by Eq. (14), and u\ is small perturbation. Sub- 
stituting this form of u along with the explicit expression 
of u in to the Eq.(16) and neglecting higher order terms 
in Mi, one would obtain the expression for u\ as follows 
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3m 
2R 2 



l + -cos2(4> + 



(17) 



For the next higher order correction term ui , one may 
follow exactly in the same way as before i.e. by writing 
u = u + u\ + «2, substituting the expressions of u and 
u\ in to Eq.(16), ignoring higher order terms in U2 and 
finally solving Eq.(16) for u 2 - The resulting solution for 
photon orbit which upto the second order accuracy in m 
is given by 



1 
R 
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(20) 



Here it is to be noted that r is smaller than R (to the 
leading order r ~ R — m) i.e. in the presence of the lens 
the distance of closest approach decreases which is quite 
natural in view of the attractive nature of the lens. But in 
that case light rays starting from the source will not reach 
the observer as the source and the observer positions are 
fixed a priori. So we have to ensure that light path is 
through both the source and the observer. When such an 
aspect is taken into consideration the parameter R in Eq. 
(18) will have to be larger than R of the orbit equation 
for undeficcted light rays (Eq. (14)). Accordingly in 
presence of the lens the distance of closest approach will 
become larger for the light trajectory between the source 
and the observer in compare to that of the undeflected 
light orbit (in the absence of the lens). 

The Eq.(18) has to be a solution of the parent Eq.(ll) 
with B = 1 — 2m i V, A = 1/B, D = 1. This requirement 
implies b = R(l + m/R) to the leading order in m. The 
Eqs.(18) and (20) completely describe the photon orbit. 
In the small angle approximation the photon orbit can be 
expressed in terms of the physically meaningful variable 
r n as 



l/r = U =^±i + ^ + ^(15^-16) + 0(m 3 ) (21) 
r ri 8ri 



B. Integration method 

Integrating the equation (11) one gets 



been expressed in terms of <j> while the reverse has been 
done in the later equation (though one may not obtain 
the one from the other easily). 



± 



A^dr 



D 2 
Bb 2 



D 



(22) 



which is also the equation of light orbit. As stated be- 
fore one has to provide the metric coefficients A, B and D 
for explicit expression of photon trajectory. However, for 
most of the relevant cases the above integration cannot 
be done analytically with any order of accuracy. Hence 
one usually expands A, B and D in terms of m/f, retain 
terms upto a desired order of accuracy and then integrate 
the above equation to obtain photon orbit [4-5]. 

To exemplify let us consider the case of the 
Schwarzschild metric. In this case the Eq.(22) is inte- 
grable and the exact solution of tfi can be described in 
terms of incomplete elliptic integral of first kind [6] but 
such a form of the solution is not very illuminating, par- 
ticularly for practical purposes. Instead usually the met- 
ric coefficients are approximated up to a certain order to 
get the approximate deflection angle. Here we are looking 
for a solution upto an accuracy of second order in m/r. 
Computing light trajectory to any given order requires 
the knowledge of every terms of the metric coefficients to 
that order in the expansion [7] . Hence in the present case 
our metric coefficients will be B — 1 — 2to /r + 2m 2 /r 2 , 
A = 1 + 2m/r + 3m 2 /(2r 2 ) and D = 1. From the fact 
that at the point of closest approach dr /d(j> must vanish, 
we get from the Eq.(ll) 



b = 



1 



TO 3 TO 2 \ ., 

^ + 2^) + °M 



(23) 



Inserting the expression of b in to the Eq. (12) and 
integrating it we get 




3to 2 / \lr 2 — r 2 1 _, 
1 :rtan 



3m 2 






tan 



\A 2 - r 2 



Vr 2 -rg(2r + r Q ) 
2r 2 (r + r ) 2 



+ O (to 3 ) 



(24) 



where £ is an integration constant. When r » r , 
the above solution becomes 



(r ->• oo) = £ 



2m 2to 



(25) 



Note that both the Eqs. (18) and (24) describe the 
same photon trajectory; in the former expression r has 



IV. DEFLECTION ANGLE 

In the Schwarzschild geometry in order not to be cap- 
tured, the apparent impact parameter b has to be at least 
greater than the critical value b c = 3v3to [4]. Here we 
are considering deflection in the weak gravity scenario. 
So in our case b is quite larger than the critical value b c . 

The geometrical configuration for the phenomenon of 
gravitational bending of light in asymptotically flat space 
time is portrayed in Figure 1. The light emitted by the 
distant source S is deviated by the gravitational source 
(Lens) L and reaches the observer O. If both the observer 
and the source are situated in the flat spacetime region 
then the angle of deflection (a = ZPQO) by the lens is 
the difference between the angle of emission (from the 
source) and the angle of reception (by the observer) with 
respect to a common polar axis. 

The difference in angular coordinates of the source 
(4>s) and the observer (<f>o) is often ascribed in the 
literature (see for instance [2-5]) as the bending angle 
{5(f) = <fis(oo) — <po{oo)) when the source and the observer 
are at infinitely large distances away from the lens. But 
how 5(f) becomes equal to the bending angle under the 
stated limit on the source/observer distance is not very 
instructive from the Figure 1. The justification for such 
equivalence is discussed below. 

In this respect the first crucial aspect is to set the polar 
axis. Note that the center of the lens has already been 
set as the origin of the co-ordinate system but so far no 
restriction has been imposed on the orientation of the po- 
lar axis. One is thus free to take any line passes through 
the point L as the polar axis. But the conventional strat- 
egy of determining bending angle does not work for an 
arbitrary orientation of the polar axis. For instance con- 
sider that the optic axis i.e. the line joining the observer 
and the lens is the polar axis (OL). Accordingly <po = 
irrespective of the distance of the observer from the lens 
whereas <f>s will be a non-zero constant unless the par- 
ticular case of perfect alignment between the source, the 
lens and the observer. Once <fis is fixed, the 5(f) i.e. the 
angular difference between the source and the observer 
would not change with the change in source distance and 
more importantly this angle certainly does not give the 
bending angle. In such a choice of the polar axis the dif- 
ference of offset angles £ at the location of the source and 
the observer should directly give the bending angle but 
there is no simple way to determine these angles. 

The standard methods require that the polar axis to be 
parallel (perpendicular) to the undcflcctcd light ray (AB 
in Fig.l). In that case the offset angle £ of Eqs. (14), 
(17), (18) , (19) becomes zero (7r/2). As r — *• oo, <f> be- 
comes very small. Denoting <fi = 4>(oo) (or 7r — 0(oo)) 
when r — > oo we get from the Eq.(18) that (f)(oo) = 



-2m I R 



157rm 



(or 7T + 2m/ R 



15-rrm \ 



The lensing 



geometry emerged from these results is shown in Fig. 2. 
Since the observer is at an infinte distance away from 
the lens, the light path LO will be parallel to OQ. Hence 
the angle 0o(oo) = ZOLA = ZOQA'. Similarly for the 
source </>s(°°) = ZSLB = ZSQB' . Hence the resulting 
deflection angle ZPQO = ZA'QP + ZB'QI = -4m/ R 
which is the standard result. The computation of bend- 
ing angle using the direct integration method is also 
based on the same principle. 

When the source and/or the observer are located at 
a finite distance away from the lens, 5<p no longer does 
represent the deflection angle. It just gives the angular 
separation between the source and the observer. This 
separation of angle is independent of the light path di- 
rection and should be the same for deflected and unde- 
tected light rays. Let us go back to the Figure 1. The 
line A'B' is parallel to AB passing through Q. So de- 
flection angle ZPQO = ZA'QO + ZB'QS. ZA'QO = 
ZKOK' — ZALO. So if tpo describes the angle that 
the tangent of the light trajectory made with at the po- 
sition of observer with the line 4>o = constant, where 
<f>o denotes the angular coordinate of the observer then 
ZA'QO = tpo- 4>o- Similarly ZB'QS = ips - <t>s- So 
the total bending angle ZPQO = i[>o + 4>s ~ 4>o — 4>s- 

The invariant formula for the cosine of the angle be- 
tween two coordinate directions P and Q 



cos(ip) 



9i j P i Q j 



(gijP'P^/HgijQ'Q^ 2 



(26) 



If P is taken as the direction of the orbit and Q is 
taken as that of the coordinate line <f> constant, then 
one may write P = (dr,dfa) — (dr/d<j>, l)d<j> (dfaltO) and 
Q = (dr,0) = (1,0) dr. Consequently for the general 
spherically symmetric space time metric given by Eq. (1) 
the angle between P and Q directions becomes [2,8] 



r\[~D 

tan(ip) = — -=-\d<j)/dr\ 



A 



The Eqs. (11) and (23) straightway give 



(27) 
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D(r)B(r ) f 2 



D(r )B(r) r 2 



-1/2 



(28) 



Substituting in to Eq.(26), we get 



tan(ip) 



D(r)B(r ) f 2 
D(r )B(f) 71 



-1/2 



(29) 



So once the metric is given, one can readily obtained 
the angle ip at any point on the light trajectory using the 
above expression. This is one of the main results of the 
present work. 



The angle that the light trajectory makes with the 
direction of the polar axis is then given by (assuming 
the gravity is very weak at the positions of the ob- 
server/source) 



ipi - fa 



(30) 



where the subscript i denotes the source, observer. The 
bending angle (a) will be then 



es + eo- 



(31) 



For e, however, the knowledge of fa are needed that 
can be obtained from the light trajectory. So for a given 
metric employing the Eqs. (26) and (27) one would obtain 
the bending angle. 

When applied to the Schwarzschild metric, we get to 
the leading order in m and l/r l 



tan(ipi) 



in 

n 



(32) 



Exploiting Eqs. (18) and (20) we finally get for small 



8r2 



2m m" 

e * = + ^T ( 157r - 16 ) 



(33) 



One may also determine the above expression using 
Eq.(24) for which one has to compute the change in an- 
gle between the source/observer and the point of closest 
approach (r ). 

Therefore the resulting deflection angle becomes 



Am 



+ —^ (157r - 16) - mr 
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d 2 

a OL 



d\s 



(34) 



where doL and dLS are the radial coordinate of the 
observer and the source respectively. The above expres- 
sion agrees with the previous results [9-13], particularly 
when doL and dz,s tend to oo. However it is worthwhile 
to note that for finite doL and dLS the above expression 
has one or two extra terms in compare to many of the 
earlier expression of the bending angle. For example, in 
[10-13] the third term of the right hand side of the above 
expression is missing which can contribute significantly 
in the solar system observations when one is looking for 
an accuracy of second order or better. Another point to 
be noted is that apparently at the second order level the 
expression given in [9] differs from the above expression. 
This is because in [9] isotropic coordinates (t, p, 6, fa) was 
used instead of the standard (Schwarzschild) coordinates 
(t, r, 9, fa) as used in this work. Hence the stated differ- 
ence can be accounted by applying the transformation 
equation between isotropic and standard coordinates, 
namely r — p(l + m/2p) 2 . 



Since r is a coordinate dependent quantity, it is not 
a measurable one and hence it is better to express the 
bending angle in terms of a proper quantity such as the 
apparent impact parameter b exploiting Eq.(23). Accord- 
ingly, the deflection angle is given by 



4m 157r m 
~b~ + ~!P r 
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d 2 
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d 2 

a LS 



(35) 



The present result contains a new term involving dis- 
tance of the source (last term of the right hand side of 
the above equation) in compare to the expression given 
in [9]. When the source will be a nearby one, such as in 
the planned astrometric missions using optical interfer- 
ometry [14], this term will also contribute significantly. 
The expression of bending angle in [9] incorporates other 
effects, such as contribution of the angular momentum 
and quadrupole moment of the lens on the bending an- 
gle, which are not considered here. 



V. DISCUSSION 

In the present work we have discussed the basic method 
of computation of gravitational bending angle and re- 
lated subtle issues to clear the prevailing confusions on 
the topic. We stressed the importance of proper choice 
of polar axis in estimating the bending angle. After de- 
veloping a general and rigorous analytical framework for 
computation of bending angle we presented a simple but 
general expression for computing bending angle without 
restricting the distances of the source and the observer 
at infinitely large level. Applying the method to the 
Schwarzschild geometry we obtained the explicit expres- 
sion of the bending angle up to an accuracy of m 2 . 

In estimating the bending angle one usually works in 
a particular coordinate system and finally the bending 



angle is presented as a function of coordinate distance 
of closest approach. For instance, if one computes the 
deflection angle in the standard coordinates, the bending 
angle is expressed in terms of r whereas working with 
isotropic coordinate the bending angle is given as a func- 
tion of p . A coordinate length like r or p is not directly 
measurable, they can only can be indirectly constructed 
from the values of actual measurements. Identifying the 
measured radius of an object with coordinate distance of 
closest approach works tolerably well only up to the first 
order level but such an approximation does not work at 
second or higher order in m [7]. So expressing bending 
angle in terms of coordinate distance of closest approach 
is not a good idea at all. One may instead express as a 
function of b since the apparent impact parameter b is a 
coordinate independent quantity. But in practical pur- 
pose b is also not suiable as it is not possible to measure 
in a real experiment. It appears that the study of bend- 
ing of a massive particle offers the only way of testing 
gravity at the second order level [7] through gravitational 
deflection measurements. 

Recently there is a serge of interest in the possibil- 
ity of measuring the cosmological constant through its 
effect on the deflection of light after Rindler and Ishak 
[8] demonstrated that cosmological constant contributes 
to the bending angle against the conventional concept of 
non-contribution. However, the expressions of bending 
angle in the presence of cosmological constant obtained 
by different authors found to differ significantly [15]. The 
method prescribed in this work has found useful appli- 
cation in this respect [16]. The present approach is also 
expected to have useful implication to the lensing phe- 
nomenon in the strong field deflection limit [17] where 
the magnitude of bending is very small. A work has been 
undertaken in this direction. 
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FIG. 1: The geometrical configuration of gravitational bend- 
ing of light. 



FIG. 2: Gravitational bending geometry when the source (S) 
and the observer (O) are at infitely large distance away from 
the lens (L). AB is the polar axis, A'B' is parallel to the polar 
axis passing through the point Q, I is the image position, R 
is a reference object 



